Inside the SU(3) quark-antiquark QCD flux tube: screening versus quantum widening 
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In lattice QCD, colour confinement manifests in flux tubes. We compute in detail the quark- 
antiquark flux tube for pure gauge SU(3) dimension D = 3 + 1 for quark-antiquark distances R 
ranging from 0.4 fm to 1.4 fm. To increase the signal over noise ratio, we apply the improved 
multihit and extended smearing techniques. We detail the gauge invariant squared components of 
the colour electric and colour magnetic fields both in the mediator plane between the static quark 
and static antiquark and in the planes of the sources. We fit the field densities with appropriate 
ansatze and we observe the screening of the colour fields in all studied planes together with the 
quantum widening of the flux tube in the mediator plane. All components squared of the colour 
fields are non-vanishing and are consistent with a penetration length A ~ 0.22 to 0.24 fm and an 
effective screening mass /j, ~ 0.8 to 0.9 GeV. The quantum widening of the flux tube is well fitted 
with a logarithmic law in R. 
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I. INTRODUCTION 

Confinement in QCD remains a central problem of 
strong interactions. It has already been established, 
both from gauge invariant lattice QCD simulations [IJ- 
3] and from experimental observations like Regge tra- 
jectories [5]-[5], that the quark-antiquark confining po- 
tential is linear, and that a flux tube develops between 
quark-antiquark static charges. Even in dynamical QCD 
where the flux tube breaks due to the creation of another 
quark and antiquark, a flux tube develops up to moder- 
ate quark-antiquark distances. Recently, the flux tubes 
have been shown to also occur in lattice QCD simula- 
tions of different exotic hadrons |10ffT4"] . Here we return 
to the fundamental quark-antiquark flux tube, to mea- 
sure in detail the profile of the SU(3) pure gauge lattice 
QCD flux tube in dimensions D = 3 + 1. We parametrize 
the flux tube profile, providing new data for a better un- 
derstanding of the confinement in QCD. 

In particular, presently two different perspectives for 
the QCD flux tube exist, possibly leading to the two dif- 
ferent flux tubes of Fig. [l] and we quantitatively compare 
them. 

Already in the 1970's, Nambu [H], *t Hooft [H] and 
Mandelstam |17| proposed that quark confinement would 
be physically interpreted using the dual version of the 
superconductivity (TBI The QCD vacuum state 

would behave like an ordinary magnetic superconductor, 
where Cooper-pair condensation leads to the Meissner ef- 
fect, and the magnetic flux is excluded or squeezed in a 
quasi-one-dimcnsional tube, the Abrikosov vortex, where 
the magnetic flux is quantized topologically. Magnetic 
charges are confined by Abrikosov-Nielsen-Olesen (ANO) 
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Figure 1: In (a) we illustrate a classical flux tube, 
similar to a solution of the Ginzburg-Landau and 
Ampere equations for a superconductor. In (b) we 
illustrate a quantum flux tube, as in a Lattice QCD 
simulation, where the widening of the flux tube occurs 
due to the zero mode string vibration. The squeezing of 
the flux tube due to the colour screening in (a) is 
masked by the widening in (b). 



vortices [20 -22] in an ordinary superconductor (Meissner 
effect). Thus, it is important for the understanding of 
confinement in QCD to measure the flux tube profile, 
and to parametrize the colour screening [23-29 . More- 
over the penetration length can be related as 

A = /i- 1 (1) 

to a possible effective mass fj, of the dual gluon, if we 
further explore the analogy between QCD and supercon- 
ductors where the field in the London equation has a 
direct relation with an effective mass of the interaction 
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Figure 2: Wilson loop and example of a plaquette for 
the computation of the electric field squared, where we 
project the D = 3 + 1 space time in a plane including 

the z axis. We compute the fields squared in three 
planes perpendicular to the charge-anticharge axis: in 
the mediator plane of the charges as illustrated in this 
Figure, in the plane of the charge and in the plane of 
the anti-charge. 



are more general that the Arvis potential, since they fit 
the D = 3 and D = 4 lattice data quite well beyond the 
tachionic distance. The Coulomb term is independent of 
the string tension a and for the physical D = 3 + 1 has 
the value — th. This is the Liischer term jJS]. The en- 
ergy spectrum of a static quark-antiquark and of its flux 
tube is certainly well defined (not tachionic) and this was 
the first evidence of flux tube vibrations found in lattice 
field theory. Moreover it was shown [45 that the width 
of the groundstate flux tube diverges when R — > oo with 
a logarithmic law, 



w 2 ~ w 2 log — (4) 



particle fields, i. e., the photon. The dual gluon mass has 
been studied by several authors, [30-37 , as well as the 
gluon effective mass, see Ref. [38] for a review of the dual 
gluon and gluon effective masses present in the literature. 
Interestingly, there is also an evidence for a gluon mass 
in the Landau Gauge [35] and in the multiplicity of par- 
ticles produced in heavy ion collisions [40]. Recently the 
penetration length started to be computed with gauge 
invariant lattice QCD techniques |38l |4"T1 14"2"] . In super- 
conductors another parameter, the coherence length £ is 
defined as well and related to the curvature of the flux 
tube profile. 

On the other hand, at quark-antiquark distances larger 
than the penetration length, the flux tube is similar to a 
quantum string. And the quantum string vibrates, even 
in the groundstate where it has zero mode vibrations. A 
fair description of the fundamental QCD flux tube - with 
charges in the triplet representation of SU(3) - is given 
by the string model, based on the Nambu-Goto Action 

EHH], 

S = -a J d 2 Z . (2) 

The energy of the quantum string with length R and 
fixed ends, with quantum transverse fluctuations quan- 
tum number n, is expressed in the Liischer term and in 
the Arvis Potential jlSlHE] . 

Vn{R) = cryj^ + ^n-^) 

7T . D - 2 N 

= CTi?+ _( n ___) + ... (3) 

In Eq. (|3j, D is the dimension of the space time. Note 
that the Arvis potential is tachionic at small distances 
since the argument of the square root is negative, more- 
over rotational invariance is only achieved for D = 26. 
Nevertheless the first two terms in the 1/R expansion 



where w 2 is the mean squared radius of the flux tube. 
This enhancement of the the flux tube transverse radius 
as R — > oo is called widening. The widening as been re- 
cently extended with two-loop calculations [17]. So far 
widening has been verified numerically for compact U(l) 
QED D = 2 + 1 lattices [48] and for non-abelian SU(2) 
D = 2 + 1 lattices @5H67]. The widening in SU(3) lat- 
tice QCD and in D — 3 + 1, which is the pure gauge 
closer to real strong interactions, has not been measured 
previously. 



In this paper, we present a SU(3) gauge independent 
lattice QCD computation in D = 3+1 for the penetration 
length and of the string quantum widening. We think 
this is a premiere both for the study of widening in SU(3) 
and for dimension as large as D = 3 + 1 . This is also the 
first attempt to separate the screening from the quantum 
widening. While the screening leads to an exponential 
decay of of the flux tube profile, the widening leads to a 
gaussian profile. 



In section II, we introduce the lattice QCD formula- 
tion. We briefly review the Wilson loop for this system, 
which was used in Bicudo et al. [68 , Cardoso et al. [69 
and Cardoso et al. |10j . and show how we compute the 
colour fields and as well as the lagrangian and energy 
densities distributions. In Section III we show the tech- 
niques we utilize to increase the signal over noise ratio. 
In Section IV we discuss our ansatz for the the width of 
the QCD flux tube. In section V, the lattice numerical 
results are shown together with their fits. Finally, we 
present the conclusion in section VI. 



3 















o 


^^^^ ^^^^ 


o 


o 


o 


o 


o 


o 


• 


o 


o o 


• 


o 


o 


o 


o 


o 


• 


o 


• o 


• 


o 


o 


o 


o 


o 


• 


o 


o o 


• 


o 


o 


o 


o 


o 


o' 


o 


•1 ° 


o 


o 


o 


o 


o 


o 












Figure 3: Left: Simple Multihit. Right: Extended 
Multihit. 
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Figure 4: Staples used in the extended spatial smearing. 



II. COMPUTATION OF THE CHROMO-FIELDS 
IN THE FLUX TUBE 

We impose our static quark-antiquark system with the 
standard Wilson W(R,T) loop [70], 



W(R,T) = Tr 



-R —T 



^(o,o, 



-R 



-1, 



f/ 4 (0,0,|,^) 



£/ 4 (0,0,|,|-l) 



^(0,0,|-l,|)...C/t( ,0,^,|) 



1/2(0,0, 



-R T 
~2~' 2 



1) 



ul(o,o, 



-R T N 
~2~' 2' 



(5) 



In the limit of large euclidean time limit T 
expectation value 



(W{R,T))=Y,\Cn\ 2 e 



oo, the 



(6) 



selects the groundstate of the static quark-antiquark sys- 
tem, aligned in the z direction with an intercharge dis- 
tance R. 

To compute the gauge invariant squared components 
of the chromoelectric and chromomagnetic fields on the 
lattice, we utilize the Wilson loop and plaquette ex- 
pectation values, 



plaquette in the lattice vertices, we average the neigh- 
bouring plaquettes. 

We define our plaquette as, 

P»» (r) = 1 - ^ReTr [U,{v)U v {v + M )£/t( r + ^(r)] , 

(8) 

which, for small lattice spacing a, can be expanded in 
a series of powers of the symmetric tensor F^ u c , which 
components are the electric and magnetic field compo- 
nents. Prior to performing the trace, the expansion reads 

P„* = 1 - ^Re Tr exp iga 2 + °^ 

L c 

= RcTrj^gV [F^F^ + 0{a)} I 

-^« 2 E[^+°( a )] TC j ( 9 ) 



(E 2 (r)) = (P(r) w ) 



(W(R,T)P(r) Ql 
(W(R,T)} 



(7) 



where the jk indices of the plaquette complement the in- 
dex i of the magnetic field. The plaquette at position 
r = (x, y, z) is computed at lattice euclidian time t = 0, 
as depicted in Fig. [2] In Eq. ^ we subtract, from the 
plaquette computed in the presence of the static charges, 
the average plaquette computed in the vacuum. This 
cancels the vacuum fluctuations of the fields. To get the 



where T c — A c /2 are the generators of the Lie algebra 
and / is the identity matrix. In abelian theories, such as 
U(l) QED, the electric and magnetic fields components 
can be computed with the plaquette at order a 2 and are 
gauge invariant. In non-abelian gauge theories, such as 
SU(3), the electric and magnetic field components are not 
gauge invariant since they depend on the colour index 
c. In SU(3) we have to go up to order a 4 to find our 
first non-vanishing gauge invariant term in the plaquette 
expansion, and it is the square of a component of the 
electric or magnetic fields. For instance E x 2 — ^2 C (E X C ) 2 
is gauge invariant, while E x c is not. Thus, to directly 
produce the squared components, we perform the trace. 

Notice the field densities defined in Eq. Q are dimen- 
sionless. To arrive at physical units , 



E 



F c F c = M 
jiv 1 fiv a 4 



l-iTr(fV) 



+ 0(a) (10) 
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Figure 5: Gap between the first excited state and the 
ground state A = V2 — V\ function of R. 



we have to multiply the dimensionless field densities by 
2/3/a 4 . 

The classical energy (H) and the lagrangian (£) den- 
sities are directly computed from the filed densities, 



(H(r)) = i((E 2 (r)) + (B 2 (r))) 
(£(r))^i((E 2 (r))-(B 2 (r))) 



(11) 
(12) 



and we can utilize any of the densities, either of the 
squared component of the fields, of the action or of the 
classical energy, to study the profiles of the flux tubes. 



Figure 6: Lagrangian density computed in the mediator 
plane for inter-charge distance of R = 12, as a function 

of the cylindrical distance r. We plot separately the 
density measured in different lines of the mediator plane 
with fixed x. At large distances r, the lattice artefacts, 
due to the square and finite lattice, produce systematic 
errors already larger than the statistical error bars 
represented in the Figure. This shows that our 
statistical noise are sufficient reduced by the extended 
multihit, the extended spatial smearing and the 
variational basis methods. 



three combined techniques, we are able to get a clear sig- 
nal, with statistical errors already smaller than the lattice 
artefacts plotted in Fig. [6j 



III. TECHNIQUES EMPLOYED TO IMPROVE 
THE SIGNAL 



To compute the static field expectation value, we plot 
the expectation value (E? (r)} or (Bf(r)) as a function of 
the temporal extent T of the Wilson loop. At sufficiently 
large T, the groundstate corresponding to the studied 
quantum numbers dominates, and the expectation value 
tends to a horizontal plateau. To compute the fields, 
we fit the horizontal plateaux obtained for each point r 
For the distances R considered, we find in the range of 
T 6 [4, 12] in lattice units, horizontal plateaux with a x 2 
/dof 6 [0.3, 2.0] . We finally compute the error bars of the 
fields with the jackknife method. 

To produce the expectation values, we utilize 1100 pure 
gauge 32 4 configurations with f3 — 6.0. This beta cor- 
responds to the lattice spacing a = 0.0983737 fm and 
a- 1 = 2.00257T GeV [73]. 

In order to reduce the noise, we utilize an improved 
version of the multihit illustrated in Fig. [3] and an ex- 
tended spatial smearing technique with staples shown in 
Fig. [1] Moreover, to reduce the contamination of the 
groundstate from excited states, we use the energy gap 
between the first excited and ground states, depicted in 
Fig. [5] calculated using a variational basis. With all 



Extended Multihit 



In the multihit [74|. 175] method we replace each tempo- 
ral link by its thermal average, with it's first neighbours 

P. A ^ * • TT ,77 / dUiUie^^^^ 

fixed, that is Ua — > Ua ~ — . 

We generalize this method by instead replacing each 
temporal link by it's thermal average with the first N 
neighbours fixed, that is, 



Ua 



By using N — 2 we are able to greatly improve the 
signal, when compared with the error reduction achieved 
with the simple multihit. Of course, this technique is 
more computer intensive than simple multihit, while be- 
ing simpler to implement than multilevel |76| and it's 
application being independent in the value of R. The 
only restriction is that R > 2N for this technique to be 
valid. 
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Figure 7: Results for minus the log of the action density in the charge mediator plane, for 
R = 4, R — 6, R — 8, R — 10, i? = 12, and i? = 14. The plots suggest the exponent is quadratic at small distances 

and linear at large distances, in agreement with our ansatz. 



B. Extended spatial smearing 

To increase the ground state overlap, we use a spatial 
extended APE, |77l 175] like smearing, namely 



SU(3) 



U 



w 2 S, 
j 



w 3 



(14) 

the staples S},, Sf, and Sf, are the ones shown in Fig. 
|4j As can be seen, this technique reduces to the common 
APE smearing when W2 = w 3 = 0. 



by the solution of the variational generalized eigensystem 
(W tf (t)>c£(t) = w n (t)(W t ,(0))c](t) (15) 

where (Wij) = (Oi(t)Oj(0)) is the correlation between 
the meson creation and annihilation operators at time t 
and in the smeared states i and j respectively. 



IV. OUR ANSATZE AND THE SEPARATION 
OF PENETRATION LENGTH, COHERENCE 
LENGTH AND QUANTUM WIDENING 

A. In the mediator plane of the two static charges 



C. Variational basis to compute A 

Even using this technique we were not able to find a 
value of t for which the plaquette to Wilson Loop cor- 
relators are stable within error bars, while still have a 
sufficiently high signal to noise ratio. To solve this, we 
note that the correlator which gives the average of field 
(F) should be given by the formula (F) t = (F) co + be~ At 
for large values of t, with A = V% — Vq, being the different 
between the first excited state which has overlap with the 
Wilson loop and the ground state potential. To compute 
A, we a use a variational basis |79l [50] of four levels of 
APE smearing, with the potentials V2 and Vq being given 



In a quantum flux tube, as in the QCD flux tube, at 
least three parameters, with the dimension of a length, 
determine the flux tube profile in the mediator plane of 
the two static charges. 

The quantum width w is a function of the flux tube 
length R and measures the widening of the flux tube due 
to the zero mode quantum oscillations of the string-like 
flux tube. 

Moreover the flux tube is not an ideal string, and it 
is due to the squeezing of the fields by the colour con- 
finement. This squeezing is expressed with two param- 
eters. The penetration length A quantifies the exponen- 
tial screening of the fields penetrating the medium. But 
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Figure 8: In (a) we the illustrate, for arbitrary parameters Eq 2 = 1, A = 1, v = 1, a = 1, our ansatz for the classical 
field and the quantum, or convoluted, one as a function of the distance to the charge axis r. In (b) the convoluted 

field and its large and small r asymptotic functions are shown. 



the flux tube cannot just be parametrized by the pene- 
tration length, because it should be differentiable at the 
centre of the flux tube, with a finite curvature. The co- 
herence length £ is related to the curvature of the field 
intensity in the centre of the flux tube. The penetra- 
tion length and the coherence length are characteristic of 
the medium (QCD in our case) where the flux tube re- 
sides, and relate to the string tension a. They should ide- 
ally be measured when the string oscillations are frozen. 
For instance, the penetration length A and the coherence 
length £ are well defined in confinement models such as 
the Ginzburg-Landau and Ampere |38| equations or in 
the Bogoliubov-de Gennes equations [55] . Notice in these 
two models for the magnetic confinement in supercon- 
ductors, the electromagnetic fields are approximated as 
classical fields, and there is no quantum widening of the 
flux tube. 

However, here we submit that with present lattice 
QCD data we can at most fit two lengths in the flux tube 
profile, because the quantum widening and the classical 
width are difficult to separate. We utilize an ansatz for 
the flux tube profile to illustrate this difficulty, and we 
work with the our convention for the cylindrical coordi- 
nates is (r, 6, z). Notice in Fig. [7] how the logarithm of 
the fields we compute is similar to a parabola at small 
distances and to a line at large distances. Thus our flux 
tube profile ansatz has an exponent inspired in the rela- 
tivistic kinetic energy, in order to interpolate between a 
gaussian at small r and an exponential decay at large r, 



F 2 (r) = F 2 exp 



v / r~ 



(16) 



iV exp 







2r 






exp 







where F corresponds to any 
nents of the squared electric or 



of the compo- 
magnetic fields 



E r 2 , E e 2 , E z 2 , B r 2 , Bg 2 , B z 2 or to the lagrangian 
density L. Our ansatz is depicted in Fig. 151 and is 
expanded for large r and for small r in Eq. ( |16[ ). Our 
ansatz is parametrized with three parameters: the 
flux tube central intensity Fb 2 , the flux tube damping 
measured by the penetration length A = and the 
flux tube central curvature radius — 2Fb 2 /(A^). In the 
Ginzburg-Landau case, the curvature is to the coherence 
length £. For a simple notation, we utilize as our 
third parameter the effective distance v. Notice the 
penetration length dominates at large distances, no 
matter how much curvature we have at the origin. 

Let us then consider a typical classical flux tube profile 
E c 2 (r) as a function of the distance r to the charge axis, 
similar to our ansatz in Eq. (16). Let us convolute the 



classical flux tube profile with a gaussian distribution, 
typical of the quantum oscillation, 



4> 2 (r) — exp 



(17) 



Notice this gaussian already has a width of w = J (r 2 ) = 
a/y/2. The result of the convolution is the quantum flux 
tube profile, 



F q u\r) 



(f> 2 (r')F cl 2 (Vr' 2 +r 2 - 2rr' cos#) dd r' dr' 



(18) 

In Fig. [8j we also show the numerical result of this inte- 
gration. 



Although we can only compute the integral in Eq. ( 18 1 
numerically for all r, we are able to compute analytically 
the profile E qu 2 {r) both close to the charge axis where 
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Figure 9: Ratios of the components of the squared fields over the Lagrangian density, for different inter-charge 
distances R. We show the ratios both computed for the field profiles computed in the mediator plane of the colour 

charges and in the planes of the colour charges. 



the profile quadratic in r, 



penetration length dominates, 



F qu \r) = F 

U- 



1 - 



A 



aX 2 + 2Xv 2 + 2a 3 2v 2 + Xa 
a\ 2 Xa 

1 



(19) 



l_«^E e (S+f) erfc(g + f)J a2 
and at large distances from the charge axis where the 



F qu 2 (r) = F 2 e Xp ( 2 ^ 2 



exp 



2r 

'T 



(20) 

These two asymptotic curves to the numerical convolu- 
tion are shown in Fig. [8] From the result of the convo- 
lution, we find that our ansatz is adequate not only for 
the fit of a classical-like flux tube, but also for the fit of 
the flux tube with quantum fluctuations, since an inter- 
polation between the two asymptotic curves yields a very 
good analytical approximation to the convolution. 
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Figure 10: Results for our fits to the profile of the action density in the mediator of the charges, for 
R = A. R = 6, R = 8, R = 10, R = 12, and R = 14. 



Importantly, the penetration length A is unaffected by 
the convolution, and it is in principle measurable at the 
long distance tail of the profile in r. However it is clear, 
both from the curvature at the origin and from the radius 
mean square, that the curvature depends on all three 
distance parameters A, v, a. Thus it is not possible, 
when error bars are significant, to separate the classical 
coherence length £ from the quantum widening a A/2. 

Moreover, with our ansatz F 2 (r) defined in Eq. 
obtain the following total width of the flux tube, consid- 
ering F 2 {r) as a distribution function, 



we 



-A 2 



\v 2 



(21) 



2 A + 2v 

Thus, our ansatz for the profile in the mediator plane is 

Table I: Fits of the profile of the flux tube, for the 
action density, in the mediator plane for the 
longitudinal component. We also consider a constant 
shift of the density, very small and not shown here 



R [a] 


10 3 £ 


A [a] 


v [a] 


X 2 /dof 


4 


3.509 ± 26.72 


2.165 ± 0.033 


0.877 ± 3.335 


4.086 


6 


2.236 ± 0.078 


2.379 ± 0.156 


2.04 ± 0.365 


2.254 


8 


1.762 ± 0.023 


2.052 ± 0.201 


4.092 ± 20.22 


1.999 


10 


1.549 ± 0.046 


2.088 ± 0.536 


5.306 ± 36.43 


1.477 


12 


1.357 ± 0.051 


0.913 ± 2.044 


17.41 ± 200.1 


1.055 


14 


1.491 ± 0.053 


0.064 ± 0.018 


268.0 ± 1392.4 


1.331 



adequate to study the total width of the flux tube as a 
function of the inter-charge distance R. In the remaining 
of this paper, we utilize Eq. [l6]to fit the profile of the flux 
tube in the mediator plane, to measure the penetration 
length A and the total widening w — -J (r 2 ). 



B. In the planes of the two static charges 

In the planes containing either the quark or the anti- 
quark static charges, only one of the three characteristic 
distances of the QCD flux tube may be measured. The 
coherence length is masked by the charges, and the quan- 
tum widening only occurs in the flux tube. Thus at most 
we may measure the screening of the Coulomb field, i e 
we can only measure the penetration length A. 

Nevertheless, for a more detailed study of the screen- 
ing, we measure the fields in planes containing one of the 
two static charges. We compare our lattice data with 
three different models for the colour fields. Without con- 
finement, one has a simple Coulomb potential, 

F qu 2 (r)=F 2 ^ , (22) 

when the distance to the charge r is smaller than the 
inter-charge distance R. If confinement does produce a 
Yukawa-like screening, the colour fields take the form, 



F qu 2 (r) = F 2 exp 



-2r 
~A~ 



Ar + 1 



(23) 
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Figure 11: Square of the width of the flux tube 
w 2 = (r 2 ) in the mediator plane, computed with our 
ansatz. The error bars arc determined with Jackknifc. 
The solid line corresponds to the fit of the widening of 
the quantum string. 
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E r ", E e \ B z z , B r and B e , are relevant to understand 
confinement. 

In Fig. [9] we show that, contrary to the dual super- 
conductor models, all components of the fields are of the 
same order of magnitude inside the flux tube. Only close 
to the charges, the larger component is E z 2 in the medi- 
ator plane and Eg 2 in the planes of the charges. 

When the distance from the charges is sufficiently 
large, all the components E 2 ~ 0.4 and all the com- 
ponents B 2 ~ —0.3 in lattice spacing units. In any case 
there is no dominant component of the colour electric 
or magnetic fields. This is an important result that any 
model of confinement should address. 

This also implies that, at sufficiently large distances 
from the charges, the parameter A and the width to, com- 
puted with any of our field densities, are essentially the 
same . 



B. Screening in the mediator plane 



Finally we may also consider a simple exponential screen- 
ing similar to the one ocurring in the mediator plane of 
the flux tube, 

F q J(r)=F 2 eiq ? (^^) , (24) 

where Fq is just a normalization parameter. 

Then it is important to check whether the penetration 
length A measured in the plane of the charges is inde- 
pendent of the charge - anticharge distance R. For a 
simple picture of the screening of the colour fields, we 
must also study if the penetration length A measured in 
the planes of the charges coincides with the penetration 
length measured in the mediator plane. 

Thus we measure the colour electric and colour mag- 
netic fields in planes including the charges. Because we 
consider long flux tubes, we choose to measure the colour 
fields in the two planes parallel to the mediator plane. 
These planes are perpendicular to the z axis, and again 
the variable measuring the distance is r — y 1 y 2 + z 2 . 



We find the noise increases with R and thus we are 
able to compute the flux tube profiles only up to R — 
14 a. We think that our noise suppression techniques are 
nevertheless sufficient, since the lattice artefacts create 
larger systematic errors than the statistical noise, see Fig. 
[6] As a word of caution we notice the systematic errors 
may contribute to increase the x 2 /dof . 

The fits of the profile of the flux tube in the mediator 
plane for the action density are shown in Fig. [10] and 
are listed in Table [T] Notice we only consider in the error 
bars the statistical error, which increases with R, thus 
decreasing the x 2 /dof with R. In the smallest distance 
R = 4 the systematic errors are larger than the statistical 
errors, and the \ 2 /dof is large. In the largest distances 
R = 12 and R = 14 the statistical errors are already 
large, and the profile parameters are not well determined. 
Nevertheless we keep this distance in our study, since the 
error in the width of the flux tube remains small up to 
R = 14. 

We remark that, although the other parameters change 
with R, the penetration length A remains the same A ~ 
2.2 a, or A ~ 0.22 fm, within the statistical error bars. 
This unique scale for the penetration length is promising 
for the theoretical understanding of confinement. 



V. FITS OF THE FLUX TUBE PROFILES 

A. The squared components of the Electric and 
Magnetic fields in both planes 

Among all densities we measure, the lagrangian or ac- 
tion density is the one with the strongest and clearest sig- 
nal, therefore this is the density we utilize to parametrize 
the profiles of the flux tube. Nevertheless all the com- 
ponents squared of the electric and magnetic fields E 2 , 



Table II: Parameters of the fits to the profile of the flux 
tube, for the action density, in the planes of the charges. 
We also consider a constant shift of the density, very 
small and not shown here. 



R[a] 


10 3 £ 


X[a] 


X 2 /dof 


4 


5.3917 ± 17.468 


2.1088 ± 0.1212 


4.8315 


6 


4.3832 ± 20.748 


2.4803 ± 0.1376 


2.1892 


8 


4.2056 ± 11.041 


2.6118 ± 0.1788 


0.9665 


12 


5.6257 ± 36.337 


2.2695 ± 0.5437 


2.5743 
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Figure 12: Lattice QCD data and fits with the exponential decay ansatz for the profile of the action density in the 

planes of the charges, shown for R = 4. R = 8 and R — 12. 



The Lagrangian density in the centre of the flux tube 
and for our largest R is of the order of 1.5 x 10 -3 in 
dimensionless units. To arrive at physical dimensions we 
have to multiply this by 2/3/a 4 = 2.5 x 10 4 GeV mi" 3 , 
and we arrive at a Lagrangian density of £ ~ 38 GeV 
fm~ 3 . 



C. Widening in the mediator plane 

Since our ansatz fits quite well the flux tube profile, 
we then utilize Eq. (21) to compute the width of the 
flux tube. Besides, we also compute the error bar or 



the width with the jacknife, method. Our results for the 
width of the flux tube in the mediator plane are shown 
in Fig. [TT] As can be seen the tube flux becomes wider 
as the quark-antiquark distance is increased. We then 
fit the flux tube width with the leading order one-loop 
computation in effective string theory |47| . corresponding 
to the linear fit, 



A + B log R 



(25) 



The fit results in, A = 0.1477 ± 0.0035 fm 2 and B = 
0.0762 ± 0.0090 fm 2 with error bars computed with jack- 
knife. Notice the error bars of the fit of the widening, for 
our larger R, are much smaller that the error bars of the 
parameters A and v of our ansatz. Nevertheless we find 
a rather small x 2 /dof — 0.383. 

The B parameter can be compared with the theoretical 
leading order [37] value for the factor of the logarithmic 
term, 



B = — — - = 0.0640028/m 2 
2tt<j 



(26) 



obtained using a string tension of y/a = 0.44 GeV |73j . 

In what concerns the constant A parameter, since it 
is positive, it is possibly larger than the corresponding 



constant of the leading order expansion of the string 
theory. Possibly this happens since the QCD flux tube 
is not tachyonic and it's width is always real and pos- 
itive. Notice a simple exponential profile, according 
to Eq. (21 1, already leads for very small distances to 
w 2 = 3)^/2 ~ 0.07/m 2 . Indeed this is similar to the 
width we get at our smaller distance of R = 4a ~ 0.4 fm. 

To comply exactly with the quantum widening of an 
infinitely thin string, the string should be much thinner 
than longer, and also much thinner than the width of the 
quantum vibrations. Indeed we have R >> A, however 
w ~ A. That our fitted factor to the logarithm is close 
to one standard deviation from the theoretical 1-loop re- 
sult, considering a large part or the width is due to the 
penetration length, is already a very interesting result. 



D. Screening in the planes of the two static charges 



We find that only one of the three ansatze in Eqs. 



(22 1, ( 23 ) and ( 24 ) fits correctly the action density in the 
planes of the charges. Both the Coulomb and Yukawa 
fields produce very poor fits of our lattice data for the 
fields. A poor fit by the Coulomb ansatz was expected 
since a flux tube is consistent with colour screening. How- 
ever the Yukawa ansatz also leads to a poor fit, and this 
indicates that the screening occurring in confinement dif- 
fers from a Yukawa screening. 

Importantly, the exponential ansatz fits correctly the 
tail of the fields in the planes of the charges, see Fig. 
[T2| Thus we have screening, though it is not a Yukawa 
screening. Moreover the fit results in a parameter A ~ 
0.22 to 0.24 fm, as listed in Table |H] . The A fitted in the 
planes of the charges is consistent with the A obtained in 
the mediator plane to the charges. 
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VI. CONCLUSIONS 

We compute the quark-antiquark flux tube in pure 
gauge SU(3) lattice QCD. We measure the profile of the 
electric and magnetic field densities both in the media- 
tor plane of the colour charges and in the planes of the 
charges. We utilize three complementary techniques to 
enhance the signal to noise ratio, and are able to reduce 
the statistical noise below the systematic errors of our 
lattice setup. 

We show the flux tube is due to screening of the elec- 
tric and magnetic field components, since we measure a 
penetration length A ~ 0.22 to 0.24 fm. The inverse of 
A may indicate an effective screening mass, possibly for 
the gluon or dual gluon, of /i ~ 0.8 to 0.9 GeV. Moreover 
the same screening parameter is universal in the sense 
it occurs in all components squared of the electric and 
magnetic fields E z 2 , E 2 , E e 2 , B z 2 , B 2 and B 2 , both 
in the mediator plane and in the charge's plane. 

However there are differences to the dual superconduc- 
tor models. The vector electric and magnetic fields are 
not gauge invariant, their squared components are the 
first gauge invariant function of the field components. 
Moreover, all the squared components have the same or- 



der of magnitude, and essentially similar profiles, thus 
the longitudinal colour electric field is not dominant. 

Importantly, this allows us to use the lagrangian den- 
sity, since it has the largest signal to noise ratio, to de- 
termine the width of the flux tube up to a distance of 
14 lattice spacings. We find that the width complies, al- 
most within one standard deviation, with the logarithmic 
widening obtained at leading order in the Nambu-Gotto 
effective string theory. 

Our results lead to a better understanding of the na- 
ture of the confining SU(3) flux tube. We hope this work 
will be useful for the theoretical understanding both of 
the QCD confinement and of string theory. 
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